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16. Abstract 


An algorithm is described for constructing the model of the 
Venus surface - as an expansion in spherical functions’. The re- 
lief expansion coefficients were obtained ’up to the coefficient 
S^. The -surface picture ^representation is given according to 

this expansion. _ The' surface model constructed was used for 
processing radar observat ions . The paper contains data show- 
ing that the use of the surface model allows improved agreement 
between the design and me.asured values -of radar ranges. 
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ANNOTATION 


An algorithm is described for constructing the model of the 
Venus surface as an expansion in spherical functions. The relief 
expansion coefficients were obtained up to the coefficient S^. The 

surface picture representation is given according to this expansion. 
The surface model constructed was used for processing radar observa- 
tions. The paper contains data showing that the use of the surface 
model allows improved agreement between the design and measured 
values of radar ranges. 
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CONSTRUCTION OP A MODEL OF THE VENUS SURFACE AND 
ITS USE IN PROCESSING RADAR OBSERVATIONS 

V.A. Borodin, V.A. Stepan'yants and V.A. Shishov 
Introduction 

Radar, ranging of Venus, Mars and Mercury has been successfully /4* 
conducted since the beginning of the 1960s. Due to improvement in 
equipment, the accuracy of radar measurements is continually increas- 
ing. The errors of the most accurate measurements of the delay time 
of a signal reflected from the surface of the planet is not over 2-3 
microseconds, which corresponds to a few hundred meters in distance. 

Radar measurements presently are widely used for the development of 
highly precise theories of motion of the solar system planets. They 
permit determination of the distance from the nearest section of the 
surface of the planet to a ground measuring point. To determine the 
location of the center of mass of a planet relative to the observation 
point, it is extremely important to take account of its topography. 

Thus, irregularities of the surface of Venus reach 5-10 km in altitude, 
which significantly affects the accuracy of determination of the dis- 
tance between the centers of mass of the earth and Venus. 

In the work of Ye.V. Pit'yeva [iJ], based on a hypsometric chart 
of the surface of Venus for approximation of the relief of the planet, 
an altitude interpolation grid has been constructed. Data on the re- 
lief were successfully used for reduction of the averaged values of the 
radar observations . 

It is advisable to present information on the surface relief of 
the planet in a form which is convenient for computer processing of 
radar observations. The volume of this information is- limited by the 
place assigned to it in the computer memory. It should also be con- 

*Numbers in the margin indicate pagination in the foreign text. 
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sidered that a radar signal is reflected from a quite extensive region 
of the planet (.approximately 100 km in diameter). Therefore, it is 
sufficient to know the average relief characteristics in processing 
radar measurements. 

With this situation taken into account, it is most convenient to 
present the description of the model of" the planet surface in the form 
of an expansion in spherical functions. Such an expansion was ob- 
tained for Mars in [6]. The task of approximating the surface relief 
of Venus in a similar manner is carried out in this work. 

The most accurate data on the surface relief of Venus was ob- 
tained with the aid of the American Pioneer-Venus-1 artificial Venus 
satellite, which was injected into a highly eccentric orbit with the 
following characteristics: 

Rotation period 24 hours 

Semimajor axis of orbit 39,500 km 

Eccentricity 0.84 

Distance to pericenter 150 km 

Distance to apocenter 66,600 km 

Pericenter longitude 17° 

Inclination of orbit 104°. 

The pericenter longitude and inclination of the orbit are reduced to 
the planet ographic coordinate system, a description of which is in 
Section 1. 

The equipment installed aboard the satellite permitted determina- 
tion of the distance from the spacecraft to the planet surface. As a 
result of processing altimeter measurements, a color topographic map 
of the surface of Venus was obtained [5]. 

The following conclusions were drawn from analysis of the ex- 
periment conducted: 



the distances of all points of the surface of Venus at which 
measurements were made from the center of mass of Venus is in the range 
from 6049 km to 6062 km; 

approximately 5 %. of the territory investigated is elevated 
more than 2 km above the average radius sphere (6051.5+0.1 km); 

approximately 60 % of the territory investigated is within 
+0.5 km of the mean level. 

Accuracy of presentation of the relief was 100-150 km horizontal- 
ly (linear resolution) and 200 m in altitude. The distances from the 
center of mass of Venus to different points on the surface of the 
planet can be recovered from the map to within the color gradations. 

This accuracy is 0.5 km in the 6049.5-6056 km distance range and 1 km 
in the 6056-6062 km range. 

1. Coordinate System 

The motions of Venus and other planets of the solar system are /6 
described in the geoequatorial coordinate system of the 1950.0 epoch, 

(XYZ) 50.0- 

For tie in of surface details of Venus, the aphrodit ographic co- 
ordinate ^system XgYgZg is used (see Fig. 1). In this system, the Zg 
axis is directed to the north pole of Venus. The Xg axis is rigidly 
tied to the surface of Venus/ and " it” is ~ in the zero meridian plane.’ 

The Yg axis supplements the coordinate system to the right. 

6 6 0 

Auxiliary planetoequatorial coordinate system X g Yg Zg also Is 
introduced. The Z D e axis of this system coindides with the Z„ axis of 

D D 

the X^Y^Zt, coordinate system. The X„ axis is directed to the point of 
o- n> o- & e 

intersection of the equators of the earth and Venus. The Yg axis sup- 
plements the system to the right. 

The transition matrix from the geoequatorial (XYZ)^q q system to 
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the aphroditographic x B Y gZg system 


P=P 2 .Pl 


is expressed by the product of transition matrix P~ from the X^Y^Z 


6 6 

2 irom tne Ag Yg Lg 
coordinate system to the XgYgZg coordinate system and transition matrix 
P^ from the (XYZ)^^ Q coordinate system to the Xg e Yg"'Zg 8 coordinate 
system. 


The elements of matrix P^ are determined by the equations [3] 

P„ »cosJ2 B 
p« - sin a* 
p a - o 

Pt, - -sln.a»cosLB 

P n m COS&bCOSLb 

Pn«sinL a 

Ps^sinaBSintB 

Pu * -cosfi| sin L a 

Pa -COSL ft , 

where fig=a 0 + 90 ° is the angular distance from the spring equinox of the 
earth (y) to the ascending node of the equator of Venus on the equator 
of the earth; ig= 90 °-S 0 is the relative inclination of the equatorial 
planes of Venus and the earth; a 0 =272.8°, 6 q= 67.2° are the coordinates 
of the North Pole of Venus [ 7 ] - 


n 


The matrix 


f sinucosu o 
j-cosw slnu o 
0 0 1 


V 


gives the rotation from the Xg Yg Zg coordinate system to the XgYgZg 
coordinate system by the angle 

_ - oj=2j 13.63° - 1 . 4814205° ' d 

between the point of intersection of the equators of the earth and 
Venus and the zero meridian of Venus [7]. In this equation, d is the 
number of Julian ephemeris days which have elapsed from the start of 
the 1950.0 epoch to the current moment. 


2. Representation of Planet Relief in the Form of Expansion in 
Spherical Functions 


The shape of the surface of the planet can be represented in the 
form of an expansion in spherical functions in the following manner 
[ 6 ] 


K- ». 



( 1 ) 


where R. is the distance from the center of mass of the planet to 
points of the surface with angular coordinates <{>^, A^, R Q is the mean /8 

radius of the planet, C , are coefficients of expansion, M' is 

the maximum power of the Legendre polynomials used in the expansion: 

ccx)- »csrci). (2) 

< ( 3 ) 

g<o : CO 

6 is the Kroneker symbol, 
mo 


Coeficients of expansion C , S nm should be selected so as to 
ensure sufficiently accurate representation of the distances from the 
center of mass of the planet to its surface. These values were ob- 
tained by processing the altimeter data of the Venus satellite, and 
it is used as the measurement material in determination of coefficients 

C , S . In this case, the number of measurements is considerably 
nm’ nm 

greater than the number of coefficients of expansion. Thus, we arrive 
at the statistical problem of determination of the parameters of the 
excess composition of the observations, which is reduced to determina- 
tion of the minimum of the functional [2] 

(5) 

to 

were -\Q(q , q 2 , . . ., q^) is the vector of the refined coefficients of 

ex'paTi^Ton7 _ R i meaS are the measured distances from the center of mass 
of Venus to its surface , q 1 = C- L0 , 3 2 =C ll> . q 3 =S n » q ii® C 20* q 5~ C 21 3 * ' ** 


s M=M' (H'+2) is the number of parameters refined, N is the 
number of observations, p t * I/S* is the weight of the i-th observation, 
3^ is the a priori observation accuracy, R. (Q) are calculated by Eq. 
(l)-(4) . 

A necessary condition of the minimum of functional (5) is satis- 
faction of the following relationships: 

w 

(Ki - Ri(Q)) - 0 ( 6 ) 


On consideration that function R (Q) is linear, system of Eq. (6) 
can be written in matrix form 

W*Q=B , (7) 


where 

R 0 £jx (sr-R.)n^ • 

i.*l ifi 

(S«nfi>cos(oXi). (sin-yOcosciXO. n'j.p; (sin^sin (i-AO. ( 8 ) 
ij (sin fv>cos (o-xo . nj «f£ (sin )cos ( 1*0 .... . nj.-ift'Csin -ft>sin (M-X.) . 

The solution of system of Eq. (7) 

Q=W -1 ' B 

is the vector of the desired coefficients of expansion. 


The values of the joined Legendre functions are calculated by 
means of the recurrent relationships [1] 

(2n + I)JC(x) « <n-m.I)P ii : i (xMn + m)Pr,(x). ( 9 ) 

Pr<x>2(m+I) ^^-lC' , (x)+(n-m)(n+»m-I)F{J , (x) = 0. ( 10 ) 


In this case, the explicit expressions of joined functions 
£<x)-r. f(x)«j3x. are obtained from Eq. 

( 2 ) — ( 4 ) . The values of P £ k (x) are calculated as a result of successive 


use of Eq. (9) and (10). 


3 . Numerical Results 


The algorithm for determination of the coefficients of expansion 
were in the form of a program in FORTRAN language. The calculations 
were carried out in a BESM-6 computer. 1944 measurements were 
processed. 

In processing, component sums (8) which correspond to the next 
measurement are determined successively. Accumulation of matrix ele- 
ments and vector B occurs as a result. The memory area required for 
solution of the problem does not exceed the memory assigned for load- 
ing the observation material, matrix W (dimension M*M), vector B (di- 
mension M), the working cells, program and results. 

The program developed can be used for representation of the sur- 
face of any planet in the form of expansion in spherical functions. 

For this, it is sufficient to assign the required mean radius of the 
planet R Q and incorporate the corresponding observations. 

Table 1 presents the results of determination of the coefficients 
of expansion in spherical functions of the surface relief of Venus up 
to harmonic inclusive. The two right hand columns contain accuracy 

characteristics of determination of coefficients -B(C nm ) and B(S nm ). 

The surface relief of Venus which corresponds to this expansion is 
represented graphically in Fig. 3. The contour lines are reckoned from 
a sphere of radius 6048 km, and they are 0.4 km apart. 

Table 2 contains the parameters which characterize the matching of 
the constructed model with the measurements: 

3 Is the root mean difference of the calculated (obtained by 
means of the Table 1 coefficients) and observed distances between the 
center of mass of Venus and points on its surface: 


/10 


TABLE 1 


n 

m 

• 10 * 

iwio ‘ 

( fcc ^ io 6 


I 

0 

23.1 


4.4 


I 

I 

- 26.8 

13.2 

2.1 

2.1 

2 

0 

- 32.7 


4.5 


2 

I 

16.2 

- II . 5 

2.2 

2.2 

2 

2 

- 26.3 

- 5.8 

2.5 

2.5 

3 

0 

. 71.0 


5.5 


3 

I 

64.3 

- 9.6 

2.6 

; 2.6 

3 

2 

2.2 

39.7 

2.4 

2.4 

3 

3 

- 3.6 

- 14.7 

2.2 

' 2.2 

4 

0 

30.4 


4.2 


4 

I 

5.1 ~ 

17. 1 

2.8 

2.8 

4 

2 

32.9 

II . 4 

2.7 

2.7 

4 

3 

- 7.0 

- 3.5 

2.2 

2.2 

4 

4 

0.4 

7.5 

1.4 

1.4 

5 

0 

15.8 


5.2 


5 

I 

36.2 

31.6 

3.2 

3.2 

5 

2 

1.8 

- 30.8 

2.7 

2.7 

5 

3 

8.1 

25.2 

2.2 

2.2 

5 

4 

6.8 

- 5.5 

2.3 

2.3 

5 

5 

- 0.5 

- 0.5 

0.8 

0.8 

6 

0 

- 12.4 


3.4 


6 

I 

26.1 

- 21.9 

3.1 

3.1 

6 

2 

5.9 

15.5 

2.8 

2.8 

6 

3 

26.0 

0.0 

2.2 

2.2 

6 

4 

- 7.0 

5.9 

2.3 

2.3 

6 

5 

4.1 

6.4 

2.2 

* 2.2 

6 

6 

- 2.4 

- 2.4 

0.8 

0.8 

7 

0 

15.2 


4.1 


7 

I 

32.2 

- 19.4 

3.2 

3.2 

7 

2 

19. 1 

2.3 

3.0 

3.0 

7 

3 

8.1 

- 15.5 

2.2 

2.2 

7 

4 

1.2 

0.0 

2.3 

2.3 

7 

5 

0.9 

- 1.2 

2.2 

2.2 

7 

6 

- 2.2 

- 5.4 

1.6 

1.6 

7 

7 

- 1.4 

- 1.5 

0.5 

0.5 

8 

0 

- 17.5 


2.8 


8 

I 

- 2.5 

- 2.9 

3.1 • 

3.1 

8 

2 

31.8 

- 13.0 

3.2 

3.2 

8 

3 

17.6 

5.4 

2.2 

' 2.2 

8 

4 

_ 8 .e 

6.3 

2.2 

2.2 

8 

5 

- 1.0 

- 3.7 

2.2 

2.2 

e 

6 

- 0.2 

- II .0 

2.0 

2.0 

8 

7 

3.3 

0.3 

1.2 

1.2 

8 

P 

- 0.8 

- 0.4 

0.2 

0.2 

9 

0 

- 16.3 


3.3 


9 

I 

9.7 

7.8 

2.9 

2.9 

9 

2 

- 25.8 

- 2 D . 2 

3.2 

3.2 

9 

3 

15.7 

12.2 

' 2.3 

2.3 

9 

4 

- 7.2 

5.4 

2.3 

2.3 

9 

5 

6.8 

7.8 

■ 2.2 

2.2 

9 

6 

- 1.7 

- 6.3 

2.2 

2.2 

9 

7 

- 3.5 

3.4 

1.7 

1.7 

9 

8 

2.8 

1.4 

0.6 

0.6 

9 

9 

0.0 

0.0 

0.0 

0.0 


A is 

observed value. 


the maximum deviation of the calculated value from the 


TABLE 2. 



6, RU 

A , 1UI 

c$epH<iecKafl uojejt 

1. 16 

8.00 

uosejib c Hcnojn>soBaHHeu 

K03$$XUHeHT0B SO S„ 

0.65 

3.31 


Key: a. Spherical model 

b. Model using coef- 
ficients up to 

Similar parameters for a spherical model of the surface of Venus 
are presented for comparison. 


4. Accounting for Asphericity of Surface of Venus in Processing Radar 
Observations 



jf 

4 


-'4 

i 


- "1 



4 


4 

1 

<* 

? 



In construction of a theory of motion of the planet with the use /13 - 
of radar measurements, the delay times of the signal reflected from the 
surface of Venus must be calculated. The distance from the center of 
mass of the planet to point C (Fig. 2) on its surface which reflects 
the radio signal (subradar point) must be known for this. In order to 
determine this distance, the planetographic coordinates of the sub- 
radar point at the time of reflection of the radio signal should be 
calculated and then, by using available data on the topography of the 
planet, the excess of the specified point of its surface above a 
sphere of a given radius should be determined. 


Let e^ and e^ be unit vectors which characterize the direction of 


the incident and reflected radio beams from the surface of the planet. 


e l =r l /r l’ 


©2 ^ 2 //r 2 5 
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where r^ are the venerocentric geoequatorial coordinates of the trans- 
mitter at signal emission time t, , r~ are the venerocentric geoequatorial- - Yl 

1 2 - q 


coordinates of the receiver at time of reception of the signal on the 
earth 


The location of the subradar point is determined by the intersec- 
tion of surface of the planet with a straight line from the center of 
mass of the planet in the direction 


Ive»l 


(ID 


In this case, angle ACE is considered to be approximately equal to 
angle BCE. 


In the aphroditographic coordinate system, unit vector of the di- 
rection to the subradar point e ’ is calculated from the relationship 

®n' =P '®n> 

where e n is the unit vector of the direction to the subradar point in 
the geoequatorial system calculated by Eq. (11) and P is the transition /lh 
matrix from the (XYZ)^ Q Q coordinate system to the X B Y B Z B coordinate 
system (see Section 1) at signal reflection time t. 

Aphroditographic latitude <p and longitude A of the subradar point 
are tied to the vector 

»®lr» • 

by the equations 

f -arc jin , X -artty (e^. /<„ ). 

Accounting for asphericity of the relief of Venus, on which the 
model of the surface of the planet constructed above Is based, ap- 
preciably decreases the discrepancies of the measured and calculated 
radar distances. 

These discrepancies are presented in Fig. 4 by longitudes of the 
subradar points on the equator, without allowance for the relief. The 
exces s o f the points of the surface above a sphere of radius 6051-5 km 
'is plotted on the same figure with a solid line. The discrepancies 


noticeably duplicate the relief of the planet. The root mean error 
is 1.24 km in this case. 

The discrepancies are represented with relief taken into account 
in Fig. 5* In this case, the root mean error is 0.99 km. Thus, allow- 
ance for the relief of Venus in processing the radar measurements pro- 
vides a significant improvement in matching the calculated and measured 
radar distances. 

It should be taken into account that the effect of the relief is 
reduced by the circumstance that the basic portion of the radar measure- 
ments is made in the area of inferior conjunctions of Venus with the 
earth (with maximum power of the earth-Venus-earth radio line). At 
inferior conjunctions, Venus turns one side to the earth which ap- 
proximately corresponds to longitude 320°. In connection with this, 
the major portion of the measurements are concentrated near this 
longitude, where the surface of Venus is quite level. 



Key: a. Zero meridian 

of Venus 

b. Equator of Venus 

c. Equator of earth 





